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ABSTRACT 

Swift triggered on a precursor to the main burst of GRB 061121 (z = 1.314), allowing observations to be made 
from the optical to gamma-ray bands. Many other telescopes, including Konus-Wind, XMM-Newton, ROTSE and 
the Faulkes Telescope North, also observed the burst. The gamma-ray. X-ray and UV/optical emission all showed a 
peak ~ 75 s after the trigger, although the optical and X-ray afterglow components also appear early on - before, 
or during, the main peak. Spectral evolution was seen throughout the burst, with the prompt emission showing a 
clear positive correlation between brightness and hardness. The Spectral Energy Distribution (SED) of the prompt 
emission, stretching from 1 eV up to 1 MeV, is very flat, with a peak in the flux density at ~ 1 keV. The optical- 
to-X-ray spectra at this time are better fitted by a broken, rather than single, power-law, similar to previous results 
for X-ray flares. The SED shows spectral hardening as the afterglow evolves with time. This behaviour might be 
a symptom of self-Comptonisation, although circumstellar densities similar to those found in the cores of molecular 
clouds would be required. The afterglow also decays too slowly to be accounted for by the standard models. Although 
the precursor and main emission show different spectral lags, both are consistent with the lag-luminosity correlation 
for long bursts. GRB 061121 is the instantaneously brightest long burst yet detected by Swift. Using a combination of 
Swift and Konus-Wind data, we estimate an isotropic energy of 2.8 x 10^'^ erg over 1 keV ~ 10 MeV in the GRB rest 
frame. A probable jet break is detected at '--^ 2 x 10*^ s, leading to an estimate of ~ 10^^ erg for the beaming-corrected 
gamma-ray energy. 

Subject headings: gamma-rays: bursts — X-rays: individual (GRB 061121) 
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Gamma- Ray Bursts (GRBs) are intrinsically extremely 
luminous objects, approaching values of lO^'' erg s~^ if 
the radiation is isotropic (e.g., Frail et al. 2001; Bloom 
et al. 2003). This energy is emitted over all bands in the 
electromagnetic spectrum; to understand GRBs as fully 
as possible, panchromatic observations are required over 
all time frames of the burst. 

The Swift multi-wavelength observatory (Gehrels et 
al. 2004) is designed to detect and follow-up GRBs. 
With its rapid slewing ability. Swift is able to follow 
bursts and their afterglows from less than a minute af- 
ter the initial trigger, and can often still detect them 
weeks, and sometimes months, later. On rare occa- 
sions, such as when Swift triggers on a precursor to the 
main burst, the prompt emission, as well as the after- 
glow, can be observed at X-ray and UV/optical wave- 
lengths. GRB 061121, the subject of this paper, is only 
the third GRB Swift has detected in this manner (af- 
ter GRB 050117 - Hill et al. 2006 and GRB 060124 - 
Romano et al. 2006), out of the almost 200 bursts trig- 
gered on in the first two years of the mission. Of 
these, GRB 061121 is the second well-sampled event 
(GRB 060124 was the first), and the first for which the 
UV/Optical Telescope (UVOT) was in event mode. 

In addition to the small number of precursor triggers, 
around 10% of S'wi/t bursts show detectable emission over 

1^ GRB 050820A would possibly have also been in this category, 
but Swift entered the South Atlantic Anomaly (SAA) just as a 
dramatic increase in count rate began (Cenko et al. 2006; Page et 
al. 2005a; Cummings et al. 2005; Page et al. 2005b; Chester et al. 
2005); Swift does not actively collect data during SAA passages. 
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the BAT bandpass by the time the narrow field instru- 
ments (NFIs) are on target. 

Besides the Swift observations of prompt emission, 
there have been a small number of prompt optical mea- 
surements of GRBs, thanks to the increasing number 
of robotic telescopes around the world. A variety of 
behaviours has been found, with some optical (and in- 
frared) light-curves tracking the gamma-ray emission 
(e.g., GRB 041219A - Vestrand et al. 2005; Blake 
et al. 2005), while others appear uncorrelated (e.g., 
GRB 990123 - Akerlof et al. 1999, Panaitescu & Kumar 
2007, though see also Tang & Zhang 2006; GRB 050904 
- Boer et al. 2006; GRB 060111B - Klotz et al. 2006; 
GRB 060124 - Romano et al. 2006). GRB 050820A (Ves- 
trand et al. 2006) showed a mixture of both correlated 
and uncorrelated optical flux. 

Where correlations exist between different energy 
bands, it is likely that there is a common origin for 
the components. In the uncorrelated cases, the opti- 
cal emission may be due to an external reverse shock 
(e.g.. Sari & Piran 1999; Meszaros & Rees 1999), while 
the prompt gamma-rays are caused by internal shocks. 
Cenko et al. (2006) suggest that the early optical data 
for GRB 050820A are produced by the forward shock 
passing through the band. In the case of GRB 990123, 
Panaitescu & Kumar (2007) have suggested that the 
gamma-rays arose from inverse Comptonisation, while 
the optical emission was due to synchrotron processes; 
they do not assume a specific mechanism for the energy 
dissipation, allowing for the possibility of either internal 
or reverse-external shocks. 

It is unclear whether precursors are ubiquitous fea- 
tures of GRBs, often remaining undetected because of 
a low signal-to-noise ratio or being outside the energy 
bandpass of the detector, or whether only some bursts 
exhibit them. A detailed discussion of the precursor phe- 
nomenon is beyond the scope of this paper and will be 
addressed in a future publication. 

In this paper, we report on the multi-wavelength ob- 
servations of both the prompt and afterglow emission of 
GRB 061121. f}2]dctails the observations made by Swift, 
Konus-Wmd, XMM- Newton, ROTSE^^ and the Faulkes 
Telescope North (FTN), with multi-band comparisons 
being made. In [JH we discuss the precursor, prompt and 
afterglow emission, with a summary given in [21 

Throughout the paper, the main burst (^ 60-200 s 
after the trigger) will be referred to as the prompt emis- 
sion, and the emission seen over —5 to -1-10 s as the 
precursor, where the BAT trigger time Tq = s. Er- 
rors are given at 90% confidence (e.g., Ax^ = 2.7 for 
one interesting parameter) unless otherwise stated, and 
the convention F^^t oc i'~^t^°' (with the photon spec- 
tral index, T = (3 + 1 where dN/dE oc E~^) has been 
followed. We have assumed a flat Universe, with Hubble 
constant. Hp = 70 km s~^ Mpc~^, cosmological constant, 
r^A = 0.73 and ^matter = 1-f^A- 

2. OBSERVATIONS AND ANALYSES 

Two years and one day after launch, the Burst Alert 
Telescope (BAT; Barthelmy et al. 2005) triggered on a 
precursor to GRB 061121 at 15:22:29 UT on 21st Novem- 
ber, 2006. Swift slewed immediately, resulting in the 

^® Robotic Optical Transient Search Experiment 



NFIs being on target and beginning to collect data 55 s 
(X-ray Telescope: XRT; Burrows et al. 2005a) and 62 s 
(UVOT; Roming et al. 2005) later. This enabled broad- 
band observations of the main burst event, which peaked 
'--^ 75 s after the trigger, leading to spectacular multi- 
wavelength coverage of the prompt emission. The most 
accurate Swift position for this burst was that determined 
by the UVOT: RA = 09'' 48"^ 54!55, decl = -13° 11' 42'.'4 
(J2000.0; 90% confidence radius of 0'.'6; Marshall et al. 
2006); the refined XRT position is only O'.'l from these 
coordinates (Page et al. 2006b). 

GRB 061121 was declared a 'burst of interest' by the 
Swift team (Gehrels et al. 2006a), to encourage an inten- 
sive ground- and space-based follow-up programme. In 
addition to the Swift observations, the prompt emission 
of GRB 061121 was detected by RHESSf^ (Bellm et 
al. 2006), Konus- Wind and Konus-A (Golenetskii et al. 
2006). Later afterglow observations were obtained in the 
X-ray [XMM-Newton - Schartel 2006) and radio (VLA^" 
- Chandra & Frail 2006) bands. ATCA^i and WSRT22 
also observed in the radio band between ~5.2 day and 
~6.2 day after the burst, but did not detect the after- 
glow (van der Horst et al. 2006a, b), implying it had faded 
since the VLA observation. 

Likewise, extensive optical follow-up observations were 
performed: ROTSE-IIIa (Yost et al. 2006), FTN (Me- 
landri et al. 2006), Kanata 1.5-m telescope (Uemura et 
al. 2006), the University of Miyazaki 30-cm telescope 
(Sonoda et al. 2006), MDM^^ (Halpern et al. 2006a,b; 
Halpern & Armstrong 2006a,b), P6024 (Cenko 2006), 
ART^s (Torii 2006), the CrAO^^ 2.6-m telescope (Efimov 
et al. 2006a,b) and SMARTS/ ANDICAM^^ (at infrared 
wavelengths, too; Cobb 2006) all detected the optical 
afterglow. Spectroscopic observations were performed 
at the Keck telescope about 12 minute after the trig- 
ger, finding a redshift of z = 1.314 for the optical af- 
terglow, based on absorption features (Perley & Bloom 
2006; Bloom et al. 2006). 

GRB 061121 has the highest instantaneous peak fiux 
of all the long bursts detected by Swift to date (e.g., 
Angelini et al. in prep). 

2.1. Gamma-ray Data 
2.1.1. BAT 

Temporal Analysis — After the initial precursor, the BAT 
count rate returned to close to the instrumental back- 
ground level, until To-l-60 s, at which point the much 
brighter main burst began. This is characterised by a 
series of overlapping peaks, each brighter than the previ- 
ous one, after which the gamma-ray flux decayed (from 
~To-l-75 s to ~To-l-140 s). Event data were collected un- 
til almost 1 ks after the trigger, thus covering the entire 
emission period. 

Reuven Ramaty High Energy Solar Spectroscopic Imager 

Very Large Array 
■^^ Australia Telescope Compact Array 
■^^ Westerbork Synthesis Radio Telescope 
'^^ Michigan-Dartmouth-MIT Observatory 
■^■^ Palomar 60 inch 
■^^ Automated Response Telescope 

Crimean Astrophysical Observatory 
■^^ Small and Moderate Aperture Research Telescope System/A 
Novel Double-Imaging CAMera 
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Tgo, over 15-150 keV, and incorporating both the pre- 
cursor and main emission, is 81 ± 5 s, measured from 
8.8-89.8 s after the trigger^^. Figure [1] shows the mask- 
weighted BAT hght-curve in the four standard energy 
bands [15-25, 25-50, 50-100, 100-150 keV; 64 ms bin- 
ning between 50-80 s after the trigger, with 1 s bins 
at all other times; units of count (fully illuminated 
detector)"^], with light-curves from other instruments: 
the precursor and the pulses of the main burst are de- 
tected over all gamma-ray bands, although the precursor 
is only marginal over the 100-150 keV BAT band. There 
is also a soft tail (detected below ~ 50 keV, when suffi- 
ciently coarse time bins are used) visible until about 140 s 
after the trigger (see bottom panel of Figure [T|) , corre- 
sponding to a similar feature in the X-ray light-curves. 

Spectral analysis — For the precursor, 
Tgo, pro = 7.7 ± 0.5 s (15-150 keV). A spectrum ex- 
tracted over this interval can be well modelled by a single 
power-law, with F = 1.68 ± 0.09 (xV^of = 26.2/23); no 
significant improvement was found by using the Band 
function (Band et al. 1993) or a cut-off power-law and a 
thermal model led to a slightly (x^ ^ 8) worse fit. The 
15-150 keV fluence for this time interval is 4 x 10^^ 
erg cm~^. 

Considering only the main event, 
Tgo, main — 18.2 ± 1.1 s (measured from 61.8-80.0 s 
post-trigger). Fitting a power-law to the mean 
spectrum during this time also results in a good fit 
(F = 1.40 ± 0.01; fluence = 1.31 x 10"^ erg cm-^ over 
15-150 keV; x^/dof =51.6/56 ); again, neither the Band 
function nor a cut-off power-law improves upon this. 
There is significant spectral evolution during the Tgo 
period, as shown in Figure [D at times when the count 
rate is higher, the spectrum is harder. This behaviour 
was also common in earlier bursts, as well as previous 
Swift detections (e.g. Golenetskii et al. 1983; Ford et al. 
1995; Borgonovo & Ryde 2001; Goad et al. 2007). The 
precursor shows a similar dependence of hardness ratio 
on count rate, suggesting that the emission processes 
in the precursor and the main burst are the same or 
similar. 

2.1.2. Konus-Wmd 

Temporal Analysis — Konus-Wind (Aptekar et al. 1995) 
triggered on the main episode of GRB 061121, while 
Konus-A triggered on the precursor (Golenetskii et 
al. 2006). Because of the spatial separation of 
Swift and Wind, the light travel-time between the 
spacecraft is 1.562 s: the Konus-Wind trigger time, 
To,K-w = To, BAT + 61.876 s. All Konus light-curves 
have been plotted with respect to the BAT trigger, 
corrected for the light travel-time. Figure [T] shows 
the Konus- Wind data plotted over the standard energy 
bands, with 64 ms binning; the bottom panel plots the 
coarser time resolution (2.944 s) 'waiting mode' data, 
showing that Konus-Wind did see slightly enhanced 
emission at the time of the precursor. The background 
levels (which have been subtracted in each case) were 
1005, 370 and 193.4 count s^^ for bands 21-83, 83-360 
and 360-1360 keV, respectively. 

■^^ Errors on the BAT Tgo are estimated to be typically 5-10%, 
depending on the shape of the light-curve. 



Spectral analysis — Table [T] gives the spectral fits to 
the Konus- Wind data in three separate time intervals 
shown by vertical lines in Figure [1] [Konus-Wind spec- 
tral intervals are automatically selected onboard): up 
to the end of the 'bump' around 70 s (the 'start' 
of the burst), the burst maximum and, finally, until 
most of the emission has died away (the burst tail). 
The data were fitted with a cut-off power-law, where 
dN/dE X e["(2-r)£;/£;p,„fc]^ leading to the photon 

indices and peak energies given in the table. The Band 
function was used to estimate upper limits for the pho- 
ton index above the peak; the values for the peak energy 
and F obtained from the Band function were the same 
as when fitting the cut-off power-law. Little variation in 
the spectral slope for energies below the peak is seen over 
these intervals, though the peak itself may have moved to 
somewhat higher energies during the burst emission. Ex- 
tracting BAT spectra over the same time intervals, and 
fitting with the same model (fixing Epoak at the value 
determined from the Konus- Wind data) results in con- 
sistent spectral indices. 

2.2. X-ray Data 
2.2.1. XRT 

Temporal Analysis — The XRT identified and centroided 
on an uncatalogued X-ray source in a 2.5 s Image Mode 
(IM) frame, as soon as the instrument was on target. 
This was quickly followed by a pseudo Piled-up Photo 
Diode (PuPD) mode frame. Following damage from a 
micrometcoroid impact in May 2005 (Abbey et al. 2005), 
the Photo Diode mode (Low Rate and Pilcd-up) has been 
disabled [see Hill et al. (2004) for details on the different 
XRT modes]; however, the XRT team are currently work- 
ing on a method to re- implement these science modes and 
to update the ground software to process the files. The 
pseudo PuPD point presented here is the first use of such 
data. 

Data were then collected in Windowed Timing (WT) 
mode starting at a count rate of ~ 1280 count s^^ (pilc- 
up corrected - see below); the rate rapidly increased to 
a maximum of ^ 2500 count s~^ at To -|- 75 s, mak- 
ing GRB 061121 the brightest burst yet detected by the 
XRT. Following this peak, the count-rate decreased, with 
a number of small flares superimposed on the underlying 
decay (see Figure [1}. Photon Counting (PC) mode was 
automatically selected when the count rate was below 
about 10 count s"'^. Around 1.5 ks, the XRT switched 
back into WT mode briefly, due to an enhanced back- 
ground linked to the sunlit Earth and a relatively high 
CCD temperature. 

Because of the high count rate, the early WT data were 
heavily piled-up; see Romano et al. (2006) for informa- 
tion about pile-up in this mode. To account for this, an 
extraction region was used which excluded the central 
20 pixels (diameter; 1 pixel = 2'.'36) and extended out 
to a total width of 60 pixels. Likewise, the ffist three 
orbits of PC data were piled-up, and the data were thus 
extracted using annular regions (inner exclusion diame- 
ter decreasing from 12 to 6 to 4 pixels as the afterglow 
faded; outer diameter 60 pixels). The count rate was 
then corrected for the excluded photons by a comparison 
of the Ancillary Response Files (ARFs) generated with 
and without a correction for the Point Spread Function 
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Fig. 1.— Top panels: Swift UVOT, XRT, BAT and Konus-Wind light-curves of GRB 061121; Icr error bars are shown for the UVOT 
and XRT data. Each instrument detected the peak of the main burst, with the precursor being detected over all gamma-ray energies. The 
vertical lines in the 360—1360 keV panel indicate the start and stop times for the spectra given in Table [Tl Bottom panel: The 15-50 keV 
BAT light-curve, with 10-s bins, showing a tail out to ~140 s. 
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start time (s) 


stop time (s) 


r 


Epeafe (keV) 


TBand 




61.876 


70.324 




478+^f 


<2.1 


72/75 


70.324 


75.188 


1 90+0. 05 
^•^"-'-0.06 




<2.9 


88/75 


75.188 


83.380 




621+282 


<2.3 


81/75 


61.876 


83.380 


1 09+U.U4 


6061^,^ 


<2.7 


95/75 



TABLE 1 

Konus-Wind cut-off power-law spectral fit results. Times are given with respect to the BAT trigger. Feand is the upper 

LIMIT OBTAINED FOR THE SPECTRAL INDEX ABOVE Epe„fc WHEN FITTING WITH THE BAND FUNCTION. 
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Fig. 2. — Top panels: Light-curves, hardness ratios (HR) and the 
variation in F using a single power-law fit during the main emission. 
The BAT light-curve (top panel) is in units of count s^^ (fully il- 
luminated detector)"^, and the corresponding hardness ratio plots 
(50-150 keV)/(15-50 keV) using 1-s binning. The XRT light-curve 
shows counts over 0.3—10 keV, while the hardness ratio compares 
(1-10 keV)/(0.3-l keV) over 1-s bins. Bottom panel: BAT hard- 
ness ratio versus count rate, showing that the emission is harder 
when brighter. Data from the precursor are shown as grey circles, 
with the main burst in black. The grey line shows a fit to the data, 
of the form HR = 0.14 CR + 0.39. 



(PSF); the ratio of these files provides an estimate of the 
correction factor. Nousek et al. (2006) give more details 
on this method. Occasionally, the afterglow was partially 
positioned over the CCD columns disabled by microm- 
eteoroid damage mentioned above. In these cases, the 
data were corrected using an exposure map. 

From To + 3 X 10^ s onwards, the afterglow had faded 
sufficiently for a nearby (41'.'5 away), constant (count 
rate ~ 0.003 count s~^) source to contaminate the GRB 
region; this source is coincident with a faint object in 
the Digitized Sky Survey and is marginally detected in 
the UVOT V filter. Thus, beyond this time, the ex- 
traction region was decreased to a diameter of 30 pix- 
els, and the count rates corrected for the loss in PSF (a 
factor of ~ 1.08). The spectrum of this nearby source 
can be modelled with a single power-law of F = 1.5^^°'^, 
with Nh = (l-8ti'2) ^ 10^^ cm~^, in comparison with 




1000 10" 10= 

time since trigger (s) 

Fig. 3.— Swift-XRT light-curve of GRB 061121. The star and 
triangle show the initial Image Mode and pseudo PuPD point (see 
text for details), followed by WT mode data (black) during the 
main burst (and at the end of the first orbit) and PC mode data 
(in grey). 
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Fig. 4. — Swift flux light-curve of GRB 061121, showing the early 
X-ray data (star, triangle and crosses) and the BAT data (grey 
histogram) extrapolated into the 0.3—10 keV band pass in units 
of erg cm~-^ together with the UVOT flux density light-curve 
(light grey circles — V-band; dark grey circles — White fllter) in units 
of erg cm~^ s~^ A~^, scaled to match the XRT flux observed at 
the start of the 'plateau' phase. 



the Galactic value in this direction of 5.09 x 10^" cm ^ 
(Dickey & Lockman 1990). 

Figure [3] shows the XRT light-curve, starting with the 
IM point (see Hill et al. 2006 for details on how IM data 
are converted to a count rate) and followed by the pseudo 
PuPD mode data. The importance of these early pre-WT 
data is clear, confirming that the XRT caught the rise of 
the main burst. 

After the bright burst, the afterglow began to follow 
the 'canonical' decay, seen in many Swift bursts (Nousek 
et al. 2006; Zhang et al. 2006a). Such a decay can 
be parameterised by a series of power-law segments; in 
this case, fitting the data beyond 200 s after the trig- 
ger (= 125 s after the main peak), two breaks in the 
light-curve were identified, with the decay starting off 
very flat (a = 0.38 ± 0.08) and eventually steepening 



to a = imtoot at 
- 32 ks (Table [2]). 



2.3 ks and then a =1.53 



+0.09 
-0.04 



at 



The addition of the second break 



GRB 061121: Broadband observations 



7 



ax 


0.38 ± 0.08 


Plateau phase 




Tbrcak, 1 


2258+377 s 






0-2 


1 n7+0-04 

^•'^'-0.06 


Shallow phase 




Tbroak,2 


(3.2lg-i) X 104 s 






az 


-I c-o+0.09 
-^•'^•-'-0.04 


Steep phase 


o 



TABLE 2 

XRT POWER-LAW LIGHT-CURVE FITS FROM 200 S AFTER THE 
TRIGGER ONWARDS; TIMES ARE REFERENCED TO THE BAT 

TRIGGER. The names used in the text for the different 

EPOCHS OF THE LIGHT-CURVE ARE LISTED IN THE LAST COLUMN. 



vastly improved the fit by Ax^ — 112.4 for two degrees 
of freedom. However, we note that O'Brien et al. (2006) 
and Wilhngale et ah (2007) advocate a different descrip- 
tion of the temporal dechne; we return to this in SJSl 

Fitting the decay of the main peak (75-200 s, keeping 
To as the trigger time) with a power-law, the slope is 
very steep, with ag = 5.1 ± 0.2. However, both Zhang 
et al. (2006a) and Liang et al. (2006) have shown that 
the appropriate time origin is the start of the last pulse. 
Thus, a model of the form f(t) oc (t— to)~"° was used, 
finding to = 58 ± 1 s and a slope of ao = 2.2+q'3; this 
is a statistically significant improvement on the power- 
law fit using the precursor To (A^^ = 32 for one extra 
parameter) . 

Figure [4] plots the Swift data in terms of flux (the BAT 
data have been extrapolated into the 0.3-10 keV band, 
using the joint fits with the XRT described in ii2.4.ip and 
fiux density for UVOT. The BAT and XRT data are fully 
consistent with each other at all overlapping times. 



Spectral Analysis — The XRT data also show that strong 
spectral evolution was present throughout the period of 
the prompt emission; this is discussed in conjunction 
with the BAT data in i j2.4.1l Considering the X-ray data 
alone, there is some indication that the spectra may be 
better modelled with a broken, rather than single, power- 
law, although the break energies cannot always be well 
constrained (see Figure [5|). For each spectrum [covering 
periods of 2 s during the main pulse, followed by two 
spectra of 5 s (80-85 s and 85-90 s) where the emission 
is fainter], the low-energy slopes were tied together for 
each spectrum (i.e., the slope measured is that averaged 
over all of the spectra), as were the high-energy indices, 
and the rest-frame column density, Nh,z, was fixed at 
(9.2 ± 1.2) X 10^1 cm-2 from the best fit to the data 
from later times (see below); only the break energy and 
the normalisation were allowed to vary. When simul- 
taneously fitting all 11 spectra, x^/dof decreased from 
142/134 to 127/132. Individually, the spectral fits were 
typically improved by of between 2-5. 

The X-ray data during the GRB 051117A flares (Goad 
et al. 2007) were found to be better modelled with bro- 
ken power-laws, with the break energy moving to harder 
energies during each flare rise, and then softening again 
as the flux decayed. Likewise, Guetta et al. (2006) found 
breaks in the X-ray spectra obtained during the flares in 
GRB 050713A. The same pattern may be occurring here, 
and there is certainly an indication of spectral curvature. 

The observed flux calculated from the spectrum corre- 
sponding to the peak of the emission (74-76 s) was mea- 
sured to be 1.66 X 10"^ erg cm~^ s~^ (over 0.3-10 keV); 
the unabsorbed value was 1.77 x 10~^ erg cm~^ s^^. 




time since trigger (s) 

Fig. 5. — Fitting the X-ray data over 0.3—10 keV with a broken 
power-law (Pi =0.69lno2 and r2 =l.Glto\t for all spectra), the 



break energy seems to move through the band, towards higher 
energies when the emission is brighter. Arrows indicate upper or 
lower 90% limits. 



The PC spectra were also extracted for the various 
phases of the light-curve ('plateau', 'shallow' and 'steep' 
- defined in Table H]) ; the results of the fitting are pre- 
sented in Table [31 In each phase, the spectrum could 
be well modelled by a single power-law (no break re- 
quired), with excess absorption in the rest-frame of the 
GRB (modelled using ZTBABS and the 'Wilms' abun- 
dance in XSPEC; Wilms et al. 2000). Together with the 
WT spectrum from ~ 200-590 s after the trigger (in the 
plateau stage), the first two PC spectra (plateau and 
shallow) are fully consistent with a constant photon in- 
dex of r = 2.07 ± 0.06 and Nh,z = (9.2 ± 1.2) x lO^i 
cm^^. 

Following the second apparent break in the light-curve, 
around 3.2 x 10''' s, the spectrum hardened slightly, to a 
1.83 ± 0.11 (or 1.87 ± 0.08 using 

-2^ 



photon index of F 
Nh z = 9.2 X 10^1 cm 



2.2.2. XMM-Newton 

XMM-Newton (Jansen et al. 2001) performed a Target 
of Opportunity observation of GRB 061121 (Observation 
ID 0311792101) less than 6.5 hr after the trigger (Schartel 
2006) and collected data for - 38 ks (MOSl, M0S2; 
Turner et al. 2001) and ~ 35 ks (PN; Striider et al. 2001). 
This observation is mainly during the 'shallow' phase, 
though also covers a short timespan after the break at 
around 32 ks. 

Figure [6] plots the PN flux light-curve and hardness 
ratio during the XMM-Newton observation, showing the 
lack of spectral evolution during this time frame; a hard- 
ness ratio calculated for the Swift data was in agreement 
with this finding. The decay slope over this time (MOSl, 
M0S2, PN and joint) is consistent with the Swift results 
(a ~ 1.3; note this crosses the time of the second break 
in the decay). 

The XMM-Newton EPIC^^ spectra show clear evidence 
for excess Nh , in agreement with the Swift data. In addi- 
tion, fitting with excess Nh in the rest-frame of the GRB 
gives a significantly better fit than at z = 0, as shown in 
Figure [T] When fitting in the observer's frame there is a 
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Epoch 


time since 
trigger (s) 


r 


Nh.z 
(10^^ cm-2) 




corresponding 
a 


Plateau 


590-1560 


2.14 ± 0.12 




62.5/52 


0.38 ± 0.08 


Shallow 


4900-22245 


2.04 ± 0.10 


8.9+li 


67.5/70 


-, „7+0.04 
J^-U'_0.06 


Steep 


34550-1152750 


1.83 ± 0.11 




48.0/55 


-'-•'J •-'-0.04 


Plateau 


590-1560 


2.09 ± 0.08 


9.2 ± 1.2 (tied) 


63.5/53 


0.38 ± 0.08 


Shallow 


4900-22245 


2.05 ± 0.06 


9.2 (tied) 


67.6/71 


1 07+°-°* 


Steep 


34550-1152750 


1.87 ± 0.08 


9.2 (tied) 


48.7/56 


1 Cq + 0.09 

^•^^•^-0.04 



TABLE 3 

XRT PC SPECTRAL FITS - REST-FRAME Nh FREE AND THEN TIED BETWEEN ALL THREE SPECTRA. THE TEMPORAL DECAY SLOPES, a, 
CORRESPONDING TO EACH STAGE ARE ALSO GIVEN. THE GALACTIC ABSORBING COLUMN OF Nh = 5.09 X 10^° CM~^ WAS ALWAYS 

INCLUDED IN THE MODEL. 




Fig. 6. — XMM- Newton EPIC-PN light-curve and hardness ratio 
of GRB 061121. The horizontal line shows the hardness ratio is 
consistent with a constant value of ~ 1.46, indicating there is no 
spectral evolution during this time. 



noticeable bump in the residuals around 0.6 keV; fitting 
with Nh at z = 1.314 removes this feature. The data 
are of sufficiently high signal-to-noise that the redshift of 
the absorber can be estimated from the spectrum. Limits 
can be placed on the redshift and absorbing column, re- 
spectively, of z > 1.2 and Nh,z > 4.6 x 10^^ cm^^ 99% 
confidence, in agreement with the spectroscopic redshift 
from Bloom et al. (2006) within the statistial uncertain- 
ties. At their value of z = 1.314, the excess Nh,z from 
the EPIC-PN spectrum is (5.3 ± 0.2) x 10^^ cm'^, lower 
than the best fit to the Swift data from the simultaneous 
'shallow' decay section, but more similar to the values ob- 
tained from fitting the optical-to-X-ray Spectral Energy 
Distributions (SEDs) in i j2.4.2l In agreement with the 
simultaneous XRT PC mode data, there is no evidence 
for a break in the EPIC spectrum over this time period. 
Spectra from neither the Reflection Grating Spectrom- 
eter (den Herder et al. 2001) nor EPIC show obvious 
absorption or emission lines. 

2.2.3. Chandra 

Chandra performed a 33 ks Target of Opportunity ob- 
servation at ~ 61 day after the trigger. No source was 
detected at the position of the X-ray afterglow, with a 
3(7 upper limit of 2.5 x 10~^^ erg cm~^ s^^. 

2.3. Optical/ UV Data 
2.3.1. UVOT 
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Fig. 7. — EPIC-PN spectrum of the late-time afterglow of 
GRB 061121, with an excess absorbing column both in the rest- 
frame of the GRB and the observer's frame. The spectrum is much 
better modelled with an excess column at z = 1.314. 

The UVOT detected an optical counterpart in the ini- 
tial White filter'^" observation, starting 62 s after the 
trigger, and subsequently in all other filters (optical and 
UV). The UVOT followed the typical sequence for GRB 
observations, with the early data being collected in event 
mode, which has a frame time of 8.3 ms during this ob- 
servation.^^ Photometric measurements were obtained 

The White filter covers a broad bandpass of A ~ 1600—6500 A. 

The data have been adjusted to take into account an incor- 
rect onboard setting (between 2006-11-10 and 2006-11-22), which 
resulted in the wrong frame times being stored in the headers of 
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from the UVOT data using a circular source extraction 
region with a 5 — 6" radius, uvotmaghist was used to 
convert count rates to magnitudes and flux; no normah- 
sation between the different filters was applied. 

As in the gamma-ray and X-ray bands, the main burst 
was detected, with an increase in count rate seen between 
~ 50-75 s after the trigger (Figures [T] and S]) . However, 
although an increase in count rate is seen for the UVOT 
data, it is by a smaller factor than observed for the XRT. 
After ~ 110 s, the UVOT emission stops decaying and re- 
brightens slightly, until 140 s after the trigger, at which 
time it flattens off and then starts to fade again (Fig- 
ure 2]). The slower decay between ~ 100-200 s may be 
indicative of the contribution of an additional (afterglow) 
component beginning to dominate. 

A single UV/optical light curve was created from all 
the UVOT filters in order to get the best measurement of 
the optical temporal decay. This was done by fitting each 
filter dataset individually (between 200 and 1 x 10^ s) 
and finding the normalisation, which was then modified 
to correspond to that of the y-band light-curve. The 
decay across all the filters beyond 200 s after the trigger 
can be fitted with a single slope of auvOT = 0.68 ± 0.02; 
the individual U, B and V decay rates are consistent 
with one another. No break in the light-curve is seen out 
to ~ 100 ks. 

2.3.2. ROTSE 

ROTSE-IIIa, at the Siding Spring Observatory in Aus- 
tralia, first imaged GRB 061121 21.6 s after the trigger 
time under poor (windy) seeing conditions. A variable 
source was immediately identified, at a position coinci- 
dent with that determined by the UVOT (Yost et al. 
2006). 

The ROTSE data (unfihcrcd, but cahbrated to the R- 
band) have been included in Figure (TU] (discussed later) . 
It is noticeable that the peak around 75 s seen in the Swift 
data is not readily apparent in the ROTSE measure- 
ments. The bandpass of the UVOT White filter is more 
sensitive to photons with wavelengths of A < 4500 A^^, 
while the ROTSE bandpass is redder. This, together 
with poor seeing conditions during the observation, may 
explain why the i?^)^^^ light-curve does not clearly show 
the main emission. 

2.3.3. Faulkes Telescope North 

The FTN, at Halcakala on Maui, Hawaii, began ob- 
servations of GRB 061121 225 s after the burst trigger, 
performing a BVRi' multi-colour sequence (Melandri 
et al. 2006). i?-band photometry was performed rela- 
tive to the USNO-B 1.0 'R2' magnitudes. Magnitudes 
were then corrected for Galactic extinction using the 
dust-extinction maps by Schlegel et al. (1998), and con- 
verted to fluxes using the absolute flux calibration from 
Fukugita et al. (1995). The photometric i?-band points 
have been included in Figure [TOl 

2.4. Broadband Modelling 

2.4.1. Gamma-rays - X-rays 

Spectral Analysis — Because the BAT was in event 
mode throughout the observation of the main burst of 

the UV OT files (Marshall 2006a). 

See |http://swiftsc. gsfc.nasa.gov/docs/heasarc/caldb/swift/ 1 



GRB 061121, detailed spectroscopy could be performed. 
Unfortunately this was not the case during the prompt 
observation of GRB 060124 (Romano et al. 2006). 

Figure [2] demonstrates the spectral evolution seen in 
both the BAT and XRT during the prompt emission. 
Spectra were extracted over 2 s intervals, in an attempt 
to obtain sufficient signal to noise while not binning over 
too much of the rapid variability. The BAT data are 
hardest around 68 s and 75 s (the second of these times 
corresponding to the peak of the main emission); the 
XRT hardness peaks about 70 s, which could be a further 
indication of the softer data lagging the harder. The joint 
spectrum (Fjoint comes from a simple absorbed power- 
law fit to the simultaneous BAT and XRT data) is at its 
hardest during the brightest part of the emission. The 
joint fit also hardens around 68-70 s, between the times 
when the BAT and XRT data respectively are at their 
hardest. The onboard spectral time-bin selection pre- 
vents the Konus- Wind data from being sliced into corre- 
sponding times, so constraints have not been placed on 
the high energy cut-off, Epoak- Breaks in the XRT-BAT 
power-laws can only be poorly constrained. 

In Figure m the BAT and XRT data were converted to 
0.3-10 keV fluxes using the time-sliced power-law fits to 
the simultaneous BAT and XRT spectra. Without the 
use of such varying conversion factors, the derived BAT 
and XRT fluxes would have been inconsistent with each 
other. 

A broadband spectrum, covering 0.3 keV to 4 MeV in 
the observer's frame (XRT, BAT and Konus- Wind) for 
~ 70-75 s post trigger was fitted by the absorbed cut-off 
power-law model described in i j2.1.2l A constant factor 
of up to 10% was included between the BAT and Konus- 
Wind data, to allow for calibration uncertainties. The 
best fit (xVdof = 301/167) gives F = 1.19 ± 0.01, with 
Epeak = GTOtlr keV. Nh,z was fixed at 9.2 x lO^i cm^^ 
(from the X-ray fits in ij2.2.1[) . Allowing F to vary be- 
tween the three spectra hints at further spectral curva- 
ture, although the differences are marginal, significant at 
only the 2a level. 

The isotropic equivalent energy (calculated using 
the time-integrated flux over the full Tgo period) is 
2.8 X 10^3 erg in the 1 keV - 10 MeV band (GRB rest 
frame), meaning that GRB 061121 is consistent with the 
Amati relationship (Amati et al. 2002). See i i3.3.2l for a 
beaming-corrected gamma-ray energy limit. 

Lag Analysts — A lag analysis (e.g., Norris et al. 1996) 
between the BAT bands leads to interesting results. 
Comparing bands 50-100 keV and 15-25 keV, the precur- 
sor emission yields a spectral lag of 600 ± 100 ms, while 
the main emission has a much smaller lag of 1 ± 6 ms. 
Note that the calculation was performed using 64 ms 
binning for the precursor and 4 ms binning for the main 
burst; see Norris (2002) and Norris & Bonnell (2006) for 
more details on the procedure. This lag for the main 
emission is rather small for a typical long burst, however 
both lags are consistent with the long-burst luminosity- 
lag relationship generally seen (Norris et al. 2000). The 
short spectral lag for the main emission, and the longer 
value for the precursor are also found when comparing 
the 100-350 keV and 25-50 keV bands. 

Similarly, comparison of the hard and soft (2-10 keV 
and 0.3-2 keV) XRT bands reveals a lag of approximately 
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Fig. 8. — Autocorrelation function of the BAT and XRT data 
during the prompt emission of GRB 061121, showing that the main 
burst peak is broader at softer energies. 



2.5 s, as the emission softens through the main burst. 
The X-ray data also lag behind the gamma-ray data, 
and the optical behind the X-ray. 

Link et al. (1993) and Fenimore et al. (1995) used a 
sample of BATSE^^ (Paciesas et al. 1999) bursts to in- 
vestigate the relationship between the duration of bursts 
and the energy band considered. They found that the 
bursts, and smaller structures within the main emission, 
generally become shorter with increasing energy (see also 
Cheng et al. 1995; Norris et al. 1996; in't Zand & Fen- 
imore 1996; Piro et al. 1998). Figure [5] plots the auto- 
correlation function over various X-ray and gamma-ray 
bands, to reinforce the point that the peak is narrower 
the harder the band - over X-ray as well as gamma-ray 
energies. Comparison of the light-curves over the differ- 
ent energy bands in Figure[T]demonstrates this as well. A 
similar behaviour was also found for GRB 060124, where 
Romano et al. (2006) compared the Tgo values obtained 
for the main burst over the X-ray and gamma-ray bands. 
Fenimore et al. (1995) found that the width of the auto- 
correlation function, W oc E~^'*, where E is the energy 
at which the function was determined; the six measure- 
ments from GRB 061121 are consistent with this finding. 

2.4.2. Optical - X-rays 

Using the Swift X-ray and UV/optical data, R and 
i' band data from the Faulkes Telescope and Rc data 
from the Kanata telescope (Uemera et al. 2006), SEDs 
were produced at epochs corresponding to the peak of 
the emission (72-75 s post BAT trigger), the plateau 
stage and during the shallow decay. Fitting at the differ- 
ent epochs gives an estimation of the broadband spectral 
variation. 

For each of the UVOT lenticular filters, the tool 
UVOt2pha was used to produce spectral files compatible 
with XSPEC, and for the latter two epochs the count rate 
in each band was set to that determined from a power-law 
fit to the individual filter light curves over the time inter- 
val in question, using a — 0.68. To determine the Faulkes 
Telescope R and i' band fiux during the plateau stage, a 
power law was fitted to the complete data set (220-1229 s 
post BAT trigger for R and 467-1401 s for i') with the 
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decay index left as a free parameter. The R magnitude 
at the mid-time of the shallow stage (6058 s) was deter- 
mined from the Kanata i?-band magnitude reported at 
6797 s (Uemera et al. 2006), assuming the same decay 
index as observed in the UVOT data. An uncertainty of 
0.2 mag was assumed as the systematic uncertainty for 
the photometric calibration of the ground based data. 

At a redshift of z = 1.314, the beginning of the Lyman- 
a forest is redshifted to an observer-frame wavelength of 
- 2812 A which falls within the UVWl filter bandpass, 
the reddest of the UV filters. A correction was applied to 
the three UV filter fluxes to account for this absorption, 
based on parameters from Madau (1995) and Madau et 
al. (1996); see also Curran et al. (in prep). 

The methods used for simultaneous fitting of the SED 
components are described in detail in Schady et al. 
(2007a). The SEDs were fitted with a power-law, or 
a broken power-law, as expected from the synchrotron 
emission, and two dust and gas components, to model 
the Galactic and host galaxy photoelectric absorption 
and dust extinction. The column density and reddening 
in the first absorption system were fixed at the Galactic 
values. [The Galactic extinction along this line of sight 
is E{B - y) = 0.046 (Schlegel et al. 1998).] 

The second photoelectric absorption system was set to 
the redshift of the GRB, and the neutral hydrogen col- 
umn density in the host galaxy was determined assuming 
Solar abundances. The dependence of dust extinction on 
wavelength in the GRB host galaxy was modelled us- 
ing three extinction laws, taken from observations of the 
Milky Way (MW), the Large Magellanic Cloud (LMC) 
and the Small Magellanic Cloud (SMC) and parame- 
terised by Pei (1992) and Cardelh et al. (1989). The 
greatest differences observed in these extinction laws are 
the amount of far UV extinction (which is greatest in 
the SMC and least in the MW) and the strength of the 
2175 A absorption feature (which is most prominent in 
the MW and neghgible in the SMC). 

Fitting these data together, a measurement of the spec- 
tral slope and optical and X-ray intrinsic extinctions (for 
the second two epochs) were obtained (Table HJ; the Ay 
values given in the table are in addition to the Ay = 0.151 
associated with the Milky Way itself. The slope above 
the break energy (which lies towards the low energy end 
of the X-ray bandpass for each phase) was assumed to 
be exactly 0.5 steeper than the spectral slope below the 
break (the condition required for a cooling break), since 
allowing all of the parameters to vary leads to uncon- 
strained fits. Figure [9] shows, as an example, the fit to 
the data in the plateau stage. 

A Milky Way dust extinction law provides the best 
overall fit to the data, using a broken power-law model, 
although the LMC model is equally acceptable. 

During the plateau phase, and adopting the bro- 
ken power-law model parameters given in Table [H we 
find gas-to-dust ratios of (1.6 ± 0.7), (2.6 ± 0.7) and 
(3.0 ± 0.7) xl022 cm-2 mag-^ for MW, LMC and SMC 
fits respectively. We can compare these estimates to the 
measured values for the MW of (4.93 ± 0.45) x 10^^ 
cm~^ mag~^ (Diplas & Savage 1994) and the LMC and 
SMC of (2.0 ± 0.8) and (4.4 ± 1.1) xlO^^ cm^^ mag~i , 
respectively (Koornneef 1982; Bouchet et al. 1985). The 
MW fit to GRB 061121, which is found to be marginally 
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Fig. 9.— Broken power-law fit to the UVOT, XRT and ground- 
based R and i' spectral energy distribution of GRB 061121 between 
~ 596-1566 s after the trigger (plateau phase) plotted in the ob- 
server's frame. The arrows indicate the beginning of the Lyman-a 
forest (I215A in the rest-frame) and the absorption feature in the 
MW dust extinction law (2175A), which is shown by a dotted line. 
The solid line corresponds to the LMC extinction, and the dashed 
one to the SMC extinction. 

the best model, is consistent with the LMC gas-to-dust 

ratio only, at the 90% confidence level. The ratios derived 
from the LMC and SMC fits are consistent with both the 
LMC and SMC gas-to-dust ratios. We note that all fits 
are inconsistent with the MW ratio at this confidence 
level, following the trend seen in pve-Swift bursts (e.g.. 
Starling et al. 2007 and references therein), and that if 
a metallicity below Solar were adopted, the gas-to-dust 
ratio of GRB 061121 would increase, moving it further 
towards the SMC value. 

3. DISCUSSION 

Swift triggered on a precursor to GRB 061121 leading 
to comprehensive broadband observations of the prompt 
emission, as well as the later afterglow. We discuss these 
here, together with possible mechanisms involved. 

3.1. Precursor 

Lazzati (2005) found that about 20% of BATSE bursts 
showed evidence for gamma-ray emission above the back- 
ground between 10 to ~200 s before the main burst, typi- 
cally with non-thermal spectra which tended to be softer 
than the main burst. GRB 060124 (Romano et al. 2006) 
and GRB 061121 show the same behaviour. 

Precursor models have been proposed for emission well- 
separated from, or just prior to, the main burst. Early 
emission occurring only a few seconds before the main 
burst has been explained by the fireball interacting with 
the massive progenitor star though the spcctnun of 
such emission is expected to be thermal (Ramirez-Ruiz 
et al. 2002a). Lazzati et al. (2007) investigated shocks in 
a cocoon around the main burst; their model predicts a 
non-thermal precursor as the jet breaks out of the sur- 
face of the star. A high-pressure cocoon is formed as the 



sub-relativistic jet head forces its way out of the star. As 
the head of the jet breaks through the surface, the energy 
of this cocoon is released through a nozzle and can give 
rise to a precursor (Ramirez- Ruiz et al. 2002a,b). Within 
the framework of this model, observers located at view- 
ing angles of 5° < 9 < 11° are expected to see first a 
relatively bright precursor, then a dark phase with lit- 
tle emission, followed, when the jet enters the imshocked 
phase, by a bright GRB; this is very similar to the light- 
curve observed for GRB 061121. Waxman & Meszaros 
(2003) demonstrate that both a series of thermal X-ray 
precursors (becoming progressively shorter and harder) 
and nonthermal emission can be produced by an emerg- 
ing shocked jot, although the nonthermal component is 
expected to be in the MeV range. There could also be 
an accompanying inverse Compton component, formed 
by the thermal X-rays being upscattered by the jet. 

The same type of smooth, wide-pulse, low intensity 
emission as seen in some precursors, but occurring af- 
ter the main emission is also occasionally seen (e.g., 
Hakkila & Giblin 2004; Nakamura 2000). Hakkila & Gib- 
lin (2004) discuss two examples where postcursor emis- 
sion is found to have a longer lag than expected from the 
lag-luminosity relation, smoother shape and to be softer. 
In the case of the GRB 061121 precursor, the spectrum is, 
indeed, softer than the main event, and shows a compar- 
atively smooth profile. The emission does have a longer 
lag than the main emission, but it is still consistent with 
the lag-luminosity relation. 

There are two expected effects which could lead to such 
a difference in lags for separate parts of a single burst: the 
much lower luminosity for the precursor (resulting from 
a much smaller Lorentz factor; the measured fluence of 
the precursor is about a factor of 30 smaller than the 
fluence of the main emission) is a natural explanation, 
while the precursor being emitted at a greater off-axis 
angle could also have an effect. In this second case, cjccta 
are considered to emerge at different angles with respect 
to the jet axis; not all of the solid angle of the jet will be 
'filled' uniformly. 

Such late postcursor emission is unlikely to be linked to 
the jet breakout from the stellar surface, and it may not 
be sensible to attribute apparently similar phenomena 
(in the form of pre- and postcursors) to entirely different 
processes. 

Pre/postcursor emission could be due to the decelera- 
tion of a faster front shell, resulting in slower shells catch- 
ing up and colliding with it (Fenimore & Ramirez-Ruiz 
1999; Umcda et al. 2005: note, however, that a faster 
shell would be inconsistent with the precursor having a 
smaller Lorentz factor as suggested to explain the lag 
discrepancy), or late activity of the central engine. The 
presence of flares in about 50% of Swift bursts is gener- 
ally attributed to continuing activity of the central en- 
gine (Burrows et al. 2005b; Zhang et al. 2006a) and the 
appearance of broken power-laws in the X-ray spectra of 
both flares and the prompt emission (Guetta et al. 2006; 
Goad et al. 2007) hints of a common mechanism. 

3.2. Prompt Emission 

The prompt emission mechanism for GRBs is still de- 
bated and the origin of Epcak is not fully understood 
(Meszaros et al. 1994; PiUa & Loeb 1998; Lloyd & Pet- 
rosian 2000; Zhang & Meszaros 2002; Rees & Meszaros 
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'^-0.15 
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0.17^^;- 
0.18«:f7 


1.22 

1.27 
1.27 


0.51 
0.72 
1.03 
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22/26 


Plateau 


PL 
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1.42 ± 0.51 


1.58 ±0.02 






0.62 ± 0.05 


167/59 






LMC 


1.98 ± 0.54 


1.64 ±0.03 






0.94 ±0.08 
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MW 
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1.39 ± 0.10 
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^^•■^^^-0.02 
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^ ^+0 3ti 

0.82«:« 

-1 r,r, + 0.25 

-■-■^^-0.20 


1.96 
2.01 
2.08 


0.52 ± 0.04 

0.74 ± 0.06 

-, no+0.09 

-■-■^•^-0.08 


84/58 
80/58 
79/58 
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PL 


SMC 


2.72 ± 0.49 


1.69 ± 0.02 






0.65 ± 0.04 


162/77 






LMC 


q 07 + 0.53 

-0.49 


1.75 ± 0.03 






n Q«+0-07 
U.98_o Qg 
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MW 




1.87 ±0.04 








127/77 




BKN PL 


SMC 
LMC 
MW 


^■'J^-0.67 
4 41+0.69 
^■^^-0.63 
4 78+°- ^5 
^- '"-0.65 


1 c:o+u.u:j 

^•°°-0.03 

1.62 ±0.03 
1.67 ±0.04 


1 30+"-^^ 

^•■^^-0.14 
, qc+0.16 
-^•■^^-0.17 


2.08 
2.12 
2.17 


0.50 ± 0.04 

0.72 ± 0.06 

1 02+°-" 
^■^^-0.10 


101/76 
99/76 
102/76 



" r2 is set to be equal to Pi + 0.5 in each broken power-law fit, as would be expected if the change in index were due to a cooling 
break. 

^ In the fit to the peak epoch, Ay is fixed to the average best-fit value found in the same model fits to plateau and shallow stage 
data. The Ay values arc given for the observer's frame of reference. 

TABLE 4 

Power-law (PL) and broken power-law (BKN PL) fits to the simultaneous UVOT and XRT spectra of GRB 061121, for 

THREE different DUST EXTINCTION MODELS: SMALL AND LARGE MAGELLANIC CLOUDS (SMC AND LMC) AND THE MiLKY WAY (MW). 
Pi AND r2 ARE THE PHOTON INDICES BELOW AND ABOVE THE SPECTRAL BREAK FOR THE BKN PL MODELS. THE DATA POINTS HAVE NOT 

BEEN CORRECTED FOR REDDENING. 



2005; Pe'er et al. 2005). The standard synchrotron 
model predicts fast cooling (Ghisellini et al. 2000) with 
a photon index, F, of 3/2 and (p/2)+l below and above 
the peak energy, respectively (e.g., Zhang & Meszaros 
2004). The Konus-Wind spectral index below Epoak is 
shallower than 3/2, which may suggest a slow cooling 
spectrum with p < 2 [Epeak being the cooling frequency 
and r =(p+l)/2] or additional heating. A slow-cooling 
spectrum can be retained by assuming that the magnetic 
fields behind the shock decay significantly in 10^-10^ cm, 
so that synchrotron emission happens in small scale mag- 
netic fields (Pe'er & Zhang 2006). 

The SED at the peak time (SED 2 in Figure [HI 
discussed below) has a peak flux density of around 
1 keV, below which the optical to X-ray spectral slope 
is 0.11 ± 0.09. This slope is harder than expected from 
the standard synchrotron model (which predicts an in- 
dex of 1/3). There should, however, be spectral cur- 
vature around the break, which could flatten the index 
(Lloyd & Petrosian 2000), so the data could still be con- 
sistent with the synchrotron model. An alternative to 
synchrotron emission, in the form of 'jitter' radiation is 
discussed by Medvedev (2000), though that model pre- 
dicts an even steeper index of 1 below the jitter break 
frequency. 

Figures [4] and \W\ show that all three instruments on- 
board Swift saw the prompt emission around 75 s after 
the BAT trigger. However, it is noticeable that most of 
the emission is in the gamma-ray and X-ray bands, with 
the optical showing a relatively small increase in bright- 
ness in comparison. Assuming the observed process is 
synchrotron, then the prompt emission which is detected 
by the UVOT will be the low-frequency extension of this 



in the internal shock. No reverse shock is apparent. 

3.3. Afterglow Emission 
3.3.1. Broken Power-law Decline Models 

The afterglow of GRB 061121 was observed over an 
even broader energy range (from radio to X-rays) than 
the prompt emission, with multi-colour data being ob- 
tained from ~ 100-10^ s after the trigger. The X-ray 
light-curve shows evidence for substantial curvature at 
later times (see Figure [3|) , as has been found for other 
Swift GRBs (e.g., GRBs 050315 - Vaughan et al. 2006; 
060614 - Gehrels et al. 2006b). The standard practice 
has been to fit such a decay using a series of power- 
law segments as a function of time. An alternative 
exponential-to-power-law description of the light-curve 
is given in i i3.3.2l 

Nousek et al. (2006) and Zhang et al. (2006a) have 
both discussed the canonical shape that many Swift af- 
terglows seem to follow: steep to plateau to shallow, 
with some light-curves showing a further steepening. In 
these previous works, the extrapolation of the BAT data 
into the XRT band was incorporated into the derivation 
of the steep decay at the start of the canonical light- 
curve shape. In the case of GRB 061121, the full curve 
can be seen entirely in X-rays, suggesting that the pre- 
vious extrapolations are reliable. For the afterglow of 
GRB 061121, only data after the end of the main burst 
have been modelled with power-laws. The early steep 
decline, which might be attributable to the curvature ef- 
fect (Kumar & Panaitescu 2000; Dermer 2004; Fan & 
Wei 2005), is not considered here. 

According to the model proposed in Nousek et al. 
(2006) and Zhang et al. (2006a), the plateau phase of 
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the light-curve is due to energy injection in the firebaU. 
The plateau phase of GRB 061121 is consistent with an 
injection of energy since the luminosity index, q, is nega- 
tive, which is the requirement for injection to modify the 
afterglow (Zhang et al. 2006a); the later two stages both 
have g > 1. However, as will be discussed in i i3.3.21 the 
plateau and final transition to the power-law decay are 
only visible in the X-ray data for GRB 061121; the start 
of the final decay is much earlier in the V and i?-bands 
(see Figure [TU|) . One might expect that energy injection 
would affect all the energy bands simultaneously, rather 
than just the X-rays. 

From the standard afterglow model computations (e.g., 
Zhang & Meszaros 2004) , we find that none of the closure 
relations fit the entire dataset completely: although the 
shallow phase (after the end of energy injection, between 
T -I- 2.3 ks and T -t- 32 ks) could be consistent with the 
evolution of a blast-wave which had already entered the 
slow cooling regime when deceleration started [i.e., v > 
max(^'m,i/c) where Vc is the cooling frequency and Vm 
is the synchrotron injection frequency; Sari et al. 1998; 
Chevalier & Li 2000] , the steeper part of the decay curve 
(T > 32 ks) is not consistent with any of the models. This 
lack of consistency suggests that a different approach is 
required. 

The change in decay slope between the shallow to 
steep phases (~ 32 ks) cannot be easily identified with 
a jet-break. It certainly seems unlikely that the simplest 
side-spreading jet model could be applicable, since the 
post-break decay index {a ^ 1.5) is not steep enough 
(a post-jet decay has a — p, where p is the electron in- 
dex). There is also some indication that the X-ray spec- 
tral slope hardens after the break, whereas no change 
in spectral signature is expected over a jet-break. In 
the case of a non-laterally expanding jet (Panaitescu 
& Meszaros 1999), a = (3/3/2) + 0.25 [for a homoge- 
neous circumstellar medium (CSM); Panaitescu et al. 

(2006) ], which does, indeed, fit the data after this break: 
[1.5 X (0.9 ± 0.08)] + 0.25 = 1.6 ± 0.1; the measured 
a is 1.53. Such a confined jet has been suggested as an 
explanation for the observed decay in a number of pre- 
vious bursts (e.g., GRB 990123 - Kulkarni et al. 1999; 
GRB 050525A - Blustin et al. 2006; GRB 061007 - 
Schady et al. 2007b). The UVOT data obtained around 
this time show little evidence for a break, whereas jet 
breaks should occur across all energy bands simultane- 
ously. However non-simultaneity could be explained by 
a multi-component outflow, where the X-ray emission is 
produced within a narrow jet, while the optical compo- 
nent comes from a wider jet with lower Lorentz factor 
(Panaitescu & Kumar 2004; Gates et al. 2007). There 
remains the issue, however, that a should steepen by 
0.75 over a jet break (Meszaros & Rees 1999), whereas 
the maximum observed change (within the 90% errors) 
is only Aa < 0.61, excluding Aa = 0.75 at almost 3cr; 
also, again there should be no spectral evolution across 
the break. There is, however, a probable jet break at 
later times, which will be covered in the next Section. 

Other multi-component models [see, e.g.. Gates et al. 

(2007) and references therein] also fail to explain the 
data, because of the lack of observed energy injection 
(plateau phase) in the optical data. 

Panaitescu et al. (2006a) discuss chromatic breaks in 
Swift light-curves, and postulate that these could be due 



to a change in microphysical parameters within a wind 
environment. However, this model requires the cooling 
frequency to lie between the X-ray and optical bands 
and, as will be discussed in !i3.3.2[ this does not seem to 
be the case here. 

3.3.2. Exponential-to-power-law Decline Model 

As first described by O'Brien et al. (2006), and further 
expanded by Willingale et al. (2007), GRB light-curves 
can be well modelled by one or two components com- 
prised of an early exponential rise followed by a power- 
law decay phase. Of these components, the first repre- 
sents the prompt gamma-ray emission and early X-ray 
decay. The second, when detected, dominates at later 
times, forming what we see as the afterglow. These re- 
sults show that fitting an intrinsically curved decay with 
multiple power-law segments runs the risk of incorrectly 
identifying temporal breaks (see also Sakamoto et al. in 
prep). In this Section the models of O'Brien et al. (2006) 
and Willingale et al. (2007) are applied to the multi-band 
afterglow data of GRB 061121. 

Figure[lO]brings together the BAT, XRT, UVOT, FTN 
and ROTSE data, along with further optical and ra- 
dio points taken from the GCN Circulars (Halpern et 
al. 2006a, b; Halpern & Armstrong 2006a, b; Chandra & 
Frail 2006; van der Horst et al. 2006a, b) and the upper 
limit from Chandra, to form a multi-energy decay plot. 
The data have been plotted as 'time since trigger -f 4 s' 
in order to include the precursor on a log time-scale. The 
optical points have all been corrected for extinction using 
Ay = 1.2 (a combination of the Galactic value of 0.151 
and an estimate of Ay ~ 1 for the GRB host galaxy - 
see !j2X2l). 

The contribution from the host galaxy reported by 
Malesani et al. (2006) and Cobb (2006) has been sub- 
tracted from the V- and i?-band flux values. The magni- 
tude of the host in the F-band is 22.4, which only changes 
the last two or three l^-band points by a small amount. 
For the i?-band we have no direct measurement, but the 
last group of MDM exposures gave an R magnitude of 
22.7, corresponding to a flux level of 2.8 /xJy, and the flux 
level is still declining at that epoch 3.3 x 10^ s), so 
an i?-band flux level of 2.5 /iJy was adopted for the host. 
The error bars shown on the last few points reflect the 
large uncertainty in the galaxy contribution subtracted. 

The curved dotted lines in Figure \W\ are the fits to 
the data using the exponential-to-power-law model, fol- 
lowed by a break to a steeper decay around 10^ s. These 
models are parameterised by the power-law decay, a, 
and Ta, the time at which this decay is established. 
For the X-ray data, Ta,x is found to be 5250l|^[! s 
and aa,x = 1-32 ± 0.03. Fits were also performed to 
the V- and i?-band data, yielding: a^y — 0.66 ± 0.04 
(with Ta,y = 70l^g s) and aa.R = 0.84 ± 0.03 

(Ta,R = 2S0tlfo s)- 

The non-detection by Chandra almost two months af- 
ter the burst shows there must have been a further steep- 
ening in the X-ray regime, and the optical data are not 
inconsistent with this finding. Constraining the tempo- 
ral index after the late break to be a = 2 (a typical slope 
for a post-jet-break decay), break times of ^ 2.5 x 10^, 
^ 2.5 X 10* and ^ 10^ s are estimated for the X-ray, V- 
and _R-band respectively; note that the UVOT y-band 
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time since trigger + 4s 

Fig. 10. — Flux density light-curves for the gamma-ray, X-ra y, op tical and radio data obtained for GRB 061121. The vertical dotted 
lines indicate the times used for the SED plots shown in Figure 1111 while the curved dotted lines show the fit to the X-ray and optical 
data, including a late-time break, as described in the text. 




Fig, 



frequency Hz 

11. — SEDs for the four time intervals indicated in Figure lTOl 



SED 2 (the peak of the burst emission) includes the Ko nus- Wind 
data, although these have not been included in Figure 1101 The 
solid lines represent the power-law fits to the BAT, XRT and Konus 
data, while the dashed lines join the radio, optical and 1 keV points. 
Spectral evolution over time is clearly seen. 



value is particularly uncertain, given the small number 
of data points at late times. Within the uncertainties, 
these times are likely to be consistent, so the turnover 
could be achromatic, as required for a jet break. From 
Willingale et al. (2007), a jet break might be expected at 
~ 100 X Ta,x - i-e., 5.5 x 10^ s, which is in agreement 
with these fits. 

As can be seen from these numbers and the mod- 
els plotted in Figure [TOl the X-ray data clearly show 
the transition from the plateau to the power-law decay, 
whereas the start of the final decay is much earlier in the 
V- and i?-bands. The l/-band decay is also significantly 



flatter (by a ^ 0.2) than that estimated for the i?-band. 
As previously stated, the V , B and U light-curves are 
all consistent with this slow decay. There have been few 
multi-colour optical decay curves obtained for GRB af- 
terglows, and, of these, the different filters [in the case of 
GRB 061007 (Schady et al. 2007; MundeU et al. 2007) 
X-ray and gamma-ray data as well as the optical] tend 
to track each other (e.g., Guidorzi et al. 2005; Blustin 
et al. 2006; de Ugarte Postigo et al. 2007). In the case 
of GRB 061121, we find that the i?-band data are fading 
more rapidly than the V . GRB 060218, which was as- 
sociated with a supernova (e.g., Campana et al. 2006a), 
shows changes throughout the optical spectra, because of 
a combination of shock break-out and radioactive heat- 
ing of the supernova ejecta. There is a large difference 
between the decays of the blue {V , U, B) and red (i?) 
data for GRB 061121, which cannot be easily explained 
by a synchrotron spectrum. Although no supernova has 
been detected in this case, we speculate that some form 
of pre-supernova thermal emission could possibly be af- 
fecting the optical data, adding energy into the blue end 
of the spectrum, thus slowing its decline. 

After the break in the decays around 10^ s, the light- 
curves across all bands become more consistent with one 
another, although there are only limited data at such a 
late time. 

The vertical dotted lines in Figure [TO] show the times 
of the SEDs plotted in Figure [TTl again, all points were 
corrected for an extinction of Ay = 1.2, so that they 
represent the true SEDs (with the frequency in the ob- 
server's frame). The solid lines represent actual fits to 
the X-ray and gamma-ray data, while the dashed lines 
just join the separate radio, optical and 1 keV points. 
The times of these SEDs, which clearly show spectral 
evolution, correspond to (1) before the main BAT peak, 
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56 s after trigger; (2) at the BAT peak, 76 s after trig- 
ger; (3) just after the start of the plateau, 300 s after the 
trigger; (4) in the main decay at 65 ks (chosen because 
radio measurements were taken at this time). SEDs 3 
and 4 do not contain any BAT or Konus data, since the 
gamma-ray flux had decayed by this point; the highest 
energy point in these corresponds to the maximum en- 
ergy (10 keV) of the X-ray fits. 

Table |4] demonstrates that the optical and X-ray spec- 
tra during the peak emission are best fitted with a broken 
power-law model, with the break energy at the very low 
energy end of the X-ray bandpass. SED 2 in Figure [TT] 
shows that this spectral break corresponds to the peak 
frequency in a flux density plot {(3i is less than zero in 
this case). Only during SED 2 is the optical flux density 
lower than that of the higher energy data. Figure [4] also 
shows that the optical emission is less strong than the 
X-ray and gamma-ray data during the main burst. 

Table [5] shows the values of a for the X-ray and opti- 
cal decays (i.e., before and after the break) in SED 4, at 
65 ks, with their corresponding spectral indices. For the 
initial stages of the power-law decay (Tq < t < 65000 s) 
the evolution of the afterglow SED and the coupling be- 
tween the temporal and spectral indices are not com- 
pletely consistent with the standard model: although the 
i?-band decay, with aa,R = 0.84 ± 0.03, is in good agree- 
ment with the homogeneous CSM model below the cool- 
ing break, the X-ray and T^-band flux decays are slower 
than expected from the measured spectral indices; they 
are in best agreement with the same constant density 
model below Uc, however. 

The point at which the power-law decay dominates 
the exponential in the optical bands is noticeably earlier 
than in the X-ray (< few hundred seconds, rather than 
~ 5000 s) and, as mentioned above, the decay indices 
are significantly different for all three (X-ray, V and R) 
bands (see Figure fTO|). At the time of SED 3, the X-ray 
data are not decaying (i.e., this is during the plateau), 
yet both the V and i?-band data have already entered 
the power-law decline phase. The i?-band is decaying 
faster than the F-band, so the spectral index through 
the optical range is becoming harder. The X-ray spec- 
tral index shows a similar hardening trend (see Table [3|), 
so the SED measured from optical to 10 keV is gradually 
getting harder. Such spectral hardening from the plateau 
to the final decay is a feature of many X-ray afterglows 
(Willingalc ct al. 2007). 

This slow hardening of the broadband spectrum with 
time could be a signature of synchrotron self-Compton 
emission (Sari & Esin 2001; Panaitescu & Kumar 2000). 
The strength of the self-Compton component in the af- 
terglow depends on the flux of low energy photons (radio- 
optical) and the electron density in the shock. Using the 
formulation in Sari & Esin (2001) the density required is 
given by 

/ fIC \ 4/3 

"1 = 3 X 10^ ^ {E,2tdayr^'^cm-^ (1) 

\ J max / 

where /ii^x//max ^he ratio of the peak flux of the 
seed synchrotron spectrum (i.e., the source of low energy 
photons) and the peak flux of the self-Compton emis- 
sion; _B52 is the isotropic burst energy in units of 10^^ erg; 
tday is the time in days after the burst (which determines 



the distance through the CSM swept up by the external 
shock). From Figure [Tl] (SEDs 1, 3 and 4) we see that 
/niax//max ^ 0.001 if the X-ray flux has a significant con- 
tribution from a self-Compton component at tday — 0.75. 
A value of £^52 = 30 gives ni k, IQ^ cm^^. Even as- 
suming the emission at 0.75 days is not dominated by 
the self-Comptonisation, and so taking the /^x//max 
ratio to be a factor of ten smaller, the density would be 
~ 5 X 10^ cm~'^, which is still high. It seems unlikely that 
self-Compton emission is the cause of the spectral hard- 
ening of the SED unless the CSM density encountered 
by the external shock is extremely large. However, there 
have been suggestions that GRBs may form in molecular 
clouds (Galama & Wijers 2001; Campana et al. 2006b, c), 
which have densities of 10** or more particles per cubic 
centimetre in the cores (Miyazaki & Tsuboi 1999; Wil- 
son et al. 1999). Typically one might expect greater red- 
dening than is found here (Table though Waxman & 
Draine (2000) discuss the possibility of dust destruction. 

The spectrum will be redshifted as the jet slows down, 
so the optical and X-ray spectral indices should, if any- 
thing, become softer - the opposite of what is seen here. 
Although spectral hardening with time is suggested from 
the data, it is not be easily explained by current models. 

Whether or not there is a Comptonised component, 
the later SEDs clearly indicate that there is a break in 
the spectrum somewhere between the optical and the X- 
ray; this is also shown by the fits in Table [H where the 
UVOT-XRT spectra are better fitted with broken power- 
laws, with Ebicak towards the low energy end of the X-ray 
bandpass. Since both the optical and X-ray bands ap- 
pear to be below the cooling frequency, from the closure 
relations given in Table [51 this change in slope cannot 
be identified with a cooling break; its origin remains un- 
clear. 

The redshift oi z — 1.314 and the isotropic energy of 
Eiso 3 X 10^^ erg ( iJ2.4.ip can be used to place con- 
straints on the jet opening angle. From Sari et al. (1999), 
and assuming that the jet break occurs at Tq -I- 2 x 10^ s, 

we have 9j ^ 4° (^)^^* ((fi)^^^ where n and 77^ are 
the density of the CSM and the efficiency of the fireball 
in converting the energy in the ejecta into gamma-rays. 
Taking 77^ = 0.2 and n = 3 cm~^ (following Ghirlanda et 
al. 2004), this gives E^ ~ 1.7 x 10^-^ erg for the beaming- 
corrected gamma-ray energy released, which is within the 
range previously determined (e.g.. Frail et al. 2001) and 
consistent with the Ghirlanda relationship (Ghirlanda et 
al. 2004). 

4. SUMMARY AND CONCLUSIONS 

Swift triggered on a precursor to GRB 061121, leading 
to unprecedented coverage of the prompt emission by all 
three instruments onboard, with the gamma-ray. X-ray 
and optical/UV bands all tracking the main peak of the 
burst. GRB 061121 is the instantaneously brightest long 
Swift burst detected thus far, both in gamma-ray and X- 
rays. The precursor and main burst show spectral lags of 
different lengths, though both are consistent with the lag- 
luminosity relation for long GRBs (Gehrels et al. 2006b). 

The SED of the prompt emission, stretching from 1 eV 
to 1 MeV shows a peak fiux density at around 1 keV and 
is harder than the standard model predicts. There is def- 
inite curvature in the spectra, with the prompt optical- 
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GRB models 


q(/3) 


"(/Sa.x)" 




"opt 

V-band ij-band 


GSM SC^ {um<u< u^) 
Wind SC= {v^<v < Uc) 
GSM or Wind SG'^ & FG'* 
(v > max(i/c, Vm)) 


■Mi 

3/3+ 1 

3/3^-1 
2 


1.49 ± 0.10 
1.99 ± 0.10 
0.99 ± 0.10 


1.32 ± 0.03 0.80 ± 0.09 
1.30 ± 0.09 
0.30 ± 0.09 


0.66 ± 0.04 0.84 ± 0.03 



Decay calculated from the measured spectral index 
^ Observed power-law decay index. 
Slow cooling. 
Fast cooling. 

TABLE 5 

Glosure relations for exponential-plus-power-law model fits to the X-ray data (/9a,x = 0.99 ± 0.07) and the 

OPTICAL-TO-X-RAY BAND (/3opt = 0.53 ± 0.06) FROM THE TIME OF SED 4 (65 KS AFTER THE BURST). 



to-X-ray spectrum being better fitted by a broken power- 
law, similar to results found for fitting X-ray flares (e.g., 
Guetta et al. 2006; Goad et al. 2007). 

The afterglow component, in both the optical and 
X-ray, starts early on ~ before, or during, the main 
burst peak (see also O'Brien et al. 2006; Willingale 
et al. 2007; Zhang et al. 2006b). The broadband 
SEDs reveal gradual spectral hardening as the afterglow 
evolves, both within the X-ray regime (F flattening from 
~ 2.05 to ~ 1.87) and between the V- and i?-band op- 
tical data (ay ~ 0.66 compared with Q!r ^ 0.84). Self- 
Comptonisation could explain the hardening, although a 
molecular-cloud-core density would be required. A prob- 
able jet-break occurs around Tq + 2 x 10^ s, shown by a 
late-time non-detection by Chandra. Before this break, 
the X-ray and y-band decays are too slow to be readily 
explained by the standard models. 

This extremely well-sampled burst shows clearly that 
there remains much work to be done in the field of GRB 
models. A single, unified model for all GRB emission 
observed should be the ultimate goal. 
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